Neutron removal cross section as a measure of neutron skin 
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We study the relation between neutron removal cross section (cr-iv) and neutron skin thickness 
for finite neutron rich nuclei using the statistical abrasion ablation (SAA) model. Different sizes 
of neutron skin are obtained by adjusting the diffuseness parameter of neutrons in the Fermi dis- 
tribution. It is demonstrated that there is a good linear correlation between ct-n and the neutron 
skin thickness for neutron rich nuclei. Further analysis suggests that the relative increase of neutron 
removal cross section could be used as a quantitative measure for the neutron skin thickness in 
neutron rich nuclei. 
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Neutron rich nuclei are expected to have a neutron skin 
due to the neutron excess and the Coulomb barrier which 
suppresses the proton density at the nuclear surface. The 
thickness of neutron skin is defined as the difference be- 
tween the root-mean-square (RMS) neutron and proton 

radii, i.e. S'„ = < >^^^ — < fp >^^^. The neutron skin 
thickness dependents on the properties of the equation 
of state (EOS) of neutron rich matter. More neutrons 
are pushed to the nuclear surface with higher pressure, 
which will form thicker neutron skin. 

It has been show that Sn is found to be related to a con- 
straint on the nuclear EOS within framework of Skyrme 
Hartree-Fock (SHF) and relativistic mean-field (RMF) 
theories [H-l^ . Strong linear correlations between Sn and 
L (the slope of symmetry energy coefficient Cgym), the 
ratio L/J (J is the symmetry energy coefficient at the 



saturation density po), J 



{O'sym is the symmetry 



energy coefficient of finite nuclei) are demonstrated 1- 
0]. The constraint on EOS is important for improving 
our knowledge of neutron rich matter and extrapolating 
the EOS to higher neutron density, and has strong im- 
pact on the property of neutron stars. It was shown that 
Sn is crucial for the transition density from nonuniform 
to uniform neutron rich matter in the crust of neutron 
star Reliable neutron densities are needed as input 
for atomic parity violation experiments Q. 

Due to the significance of Sn for studying structure of 
neutron rich nuclei, property of neutron rich matter and 
neutron star etc., measurement of Sn has become a very 
important subject in the forefront of nuclear physics re- 
search. It is remarkable that such a measurement is so 
important for investigating the property of very small 
entities like atomic nuclei as well as very big objects like 
neutron stars. As we know, the proton radius of the 
nucleus can be determined with a very high accuracy 
through the electron scattering measurement [9|. How- 
ever, the experimental accuracy of neutron radius is much 
lower compared with that of the proton radius so far. 

Recently, some methods have been proposed or used 
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to measure the radius of neutron distribution. A model 
independent method using parity violating electron scat- 
tering to measure the neutron radius of ^"^Pb to 1% ac- 
curacy (±0.05 fm) was proposed at JLab [l^. There are 
methods to extract neutron densities from proton-nucleus 
clastic scattering and antiprotonic atom data [ll-l^]- 
Another method of using the relationship between the 
nuclear matter radius and the reaction cross section {afi) 
was shown [T^. But it is also useful to have alternative 
methods of obtaining Sn directly. The dependence of 
neutron removal cross section (cr_Ar) on the neutron skin 
thickness was studied and strong linear correlation was 
observed fld\. In this work, it is demonstrated that the 
relative increase of neutron removal cross section could 
be used as a quantitative measure for Sn by using the 
statistical abrasion ablation (SAA) model . 

The SAA model was developed by Brohm et al to de- 
scribe heavy ion collisions at high energies in a p icture 
of quasi- free nucleon-nucleon (N-N) collisions [17]. This 
model takes independent N-N collisions for participants 
in an overlap zone of the two colliding nuclei and de- 
termines the distributions of abraded neutrons and pro- 
tons. The spectators outside the overlapping region are 
left to move almost without being disturbed. The proba- 
bility of abrading zero nucleons is calculated and the re- 
sulting energy-dependent reaction cross section turns out 
to be equivalent to the successful microscopic Glauber- 
type calculation. It has been demonstrated that the SAA 
model can reproduce the experimental fragment produc- 
tion cross section of heavy ion collisions very well [3 
and is useful for studying the isospin effect and isoscal- 
ing behavior of projectile-like fragments in heavy 

ions reactions. Here we give a brief description of the 
SAA model. Its details can be found in Ref. [l7| . 

In SAA model, nuclear reaction is described as two 
stages which occur in two different time scales. The 
first stage is abrasion process which describes the pro- 
duction of pre-fragment with certain excitation energy 
{E*) through the N-N collisions in the overlap zone of 
the colliding nuclei. The collisions are described by a pic- 
ture of interacting tubes. Assuming a binomial distribu- 
tion for the absorbed neutrons and protons of projectile 
in the interaction of a specific pair of tubes, the distri- 
butions of the total abraded neutrons and protons are 
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determined. For an infinitesimal tube in the projectile, 
the transmission probabilities for neutrons (protons) at 

a given impact parameter b (= |6|) are calculated by 



t]^{f-b) = exp{ 



&)^pk]}' (1) 



T 

where is thickness function of the target, which is 
normalized by / (fr'D'^{r') = and / (Pr'Dj{r') = 

Z with N and Z referring to the neutron and proton 
number in the target respectively, the vectors r, r' and 
b are defined in the plane perpendicular to the beam, 
and (y^k ^^e free nucleon-nucleon cross sections (k', 
k=n for neutron and k', k=p for proton). The thickness 
function of the target is given by 
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(2) 



with p-^ being the neutron (proton) density distribution 
of the target. The thickness function of the projectile 
D^(r') is defined in the same way. When the thickness 

function is used in the calculation, the variable r' equals 
f—b for the target and r' equals r for the projectile. So 
the average abraded mass at a given impact parameter b 
is calculated by the expression 



{^A{b)) = {^N{b)) + (Az(6)) 

= !<PrDl{r)[l-tn{r-b)] 
+ /dHi?P(f)[l-tp(f-6)], 



(3) 



where (AiV(5)) and {AZ{b)) are the average abraded 
neutron and proton number at impact parameter 6, is 
thickness function of the projectile. The integrals extend 
over the plane perpendicular to the beam direction. In 
a similar way, the variance of can be determined by 
summing the variances of all the binomial distributions 
for the tubes pairs. 



{[AA{b) - {AA{b))r) 



d^rD^{r}[l-tn{r-b)]tn{f-b) 



d\D^if)[l ~ tpif- 6)]tp(f - 6), 



(4) 



Higher moments of the distributions of abraded neutrons 
and protons could be obtained similarly. However, for the 
present purpose, sufficient information is provided by the 
first and second moments. In addition, the probability for 
zero nucleon loss is important because it is related to the 
reaction cross section. The deduced total transmission 
probability T{b) and reaction cross section are [l3| 



r(6) = exp(- / d'r{Dl{f)Djir-b)ann 



-b) 



D^if)D'^if-b)]anp 



(5) 



j dH[l-T{b)], 



(6) 



The production cross section for a specific fragment 
can be calculated from 



(j{AN,AZ) = J d^bP{AN,b)PiAZ,b), 



(7) 



where AN, AZ are the number of abraded neutrons and 
protons from the projectile. P{AN, b) and P{AZ, b) are 
the probability distributions for the abraded neutrons 
and protons at a given impact parameter 6, respectively. 
These probability distributions are determined by a su- 
perposition of different binomial distributions. 

The second stage is evaporation process in which the 
system reorganizes due to excitation. It deexcites and 
thermalizes by emission of neutron, proton and other 
light particles using the statistical model [1^. The ex- 
citation energy for projectile spectator is estimated by a 
simple relation of E* = 13.3(A^(fe)) MeV where 13.3 is 
a mean E* for an abraded nucleon from the projectile. 
This E* was given by the statistical hole-energy model 
as described in Ref. (2^]. After the evaporation stage, we 
can obtain the final fragments which are comparable to 
the experimental data. 

For the density distributions of the projectile and tar- 
get, harmonic-oscillator (HO) distributions are suitable 
for light nuclei but it is more reasonable to use the Fermi 
distribution for heavier mass system (23j . In our calcu- 
lations, the Fermi distributions are used for neutron and 
proton, respectively. They have the following form 
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-cxp(t-7o 



n,P, 



(8) 



where is the normalization constant which ensures 
that the integration of the density distribution equals 
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FIG. 1: (Color online). Correlation between cr-jv and Sn for 
Na, P, Ca and Ni isotopes with different N ~ Z on^^C target 
at lA GeV calculated by the SAA model. The straight lines 
just guide the eye. 
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the numbers of neutron («=n) and proton (i=p); U is the 
diffuseness parameter; Q is the half the density radii and 
determined by droplet model [13] 



a^N with respect to cr° ^, we define a difference factor 



a^R,[l-{h/R,)%i = n,p, 



(9) 



with bi = 0.413ii, Ri being the equivalent sharp surface 
radii of neutron and proton. Ri and ti are given by the 
droplet model. 

Neutron removal cross section (tT_7v) is defined as the 
summation of production cross sections for fragments 
with the same proton number as the projectile, but less 
neutron number than the projectile, i.e. collision pro- 
cesses with one or more neutrons being removed from 
the projectile: 



N 

= ^a(z,0), 

1=1 



(10) 



where N is the neutron number of the projectile. a{i, 0) 
is the fragment production cross section with AA^ — i 
and AZ — given by Eq. ([7]). It should be pointed out 
that (t(«,0) decreases very quickly as i increases. The 
value of cr(iV, 0) will be very close to zero for most pro- 
jectiles. Since the charge number of the fragments equals 
to that of the projectile, it indicates that peripheral reac- 
tions will be the main contributions of ct-at. Thus cr_jv 
maybe sensitive to the neutron distribution in the nu- 
clear surface. Since the reaction cross section reflects the 
information about the overall size of a nucleus, Sn of Na 
isotope determined by the interaction cross section (ct/) 
has an error in the order of ±0.06 ^ ±0.3fm when the 
error of cr/ is around 1^5% [15]. Below we will demon- 
strate that determining Sn through cr_jv measurement 
may reduce the uncertainty. 

In Ref. [l3|, Trzcihska et al found that the half den- 
sity radii for neutrons and protons in heavy nuclei were 
almost the same, but the diffuseness parameter for neu- 
trons was larger than that for protons. To study the 
dependence between g-n and Sn, the neutron diffuse- 
ness parameter in the droplet model is changed to obtain 
different value of Sn for a nucleus. Six values of Sn vary- 
ing from 0.0 fm to 0.5 fm with step of 0.1 fm are chosen. 
Different values of a^M are obtained in the SAA calcu- 
lation with different 5„ of the projectile. Thus we can 
obtain the correlation between cr_Ar and Sn- 

In this work, calculations have been done for 

23,25,27,29,31,33]\^g^ 31, 33, 35, 37, 39, 41p 42, 44, 46, 48, 50, 52^^^ ^^^^ 

58,60,62,64,66,68Nj projcctiles On ^^'c target at lA GeV. 
The results are shown in Fig.[T] When Sn varies from 0.0 
fm to 0.5 fm, we can see that (T_jv is increasing linearly 
for all the isotopes calculated. As its definition shows, 
Sn is determined by the radii of protons and neutrons. 
Since most nuclei have precise proton radii, the value of 
neutron removal cross section for Sn = or r„ = 

jv) could be calculated by the SAA model using the 
experimental rp data . To study the relative increase of 



-N 



(11) 



The systematic behavior of d with respect to Sn for all 
the isotopes are calculated and shown in Fig. [51 From 
this figure, we found that the values of d for one iso- 
tope is almost constant as the horizontal lines showing, 
especially for 5„ < 0.4 fm. For one isotope with same 
neutron skin thickness, the mean uncertainty of d is less 
than ±1.2% for Sn < 0.4 fm. It indicates that the relative 
increases of (t_jv with respect to cr^ijv almost same 
for one isotope having different neutron excess {N — Z) 
when the neutron skin thickness changes from fm to 
0.4 fm. The mean values of d for one isotope increase as 
Sn increases. We denote this mean value of d by d. The 
dependence of d on Sn for the four isotopes are shown in 
Fig. 3. Clear linear relationship between d and Sn can be 
seen. We fitted the calculated results by using a function 
of d = slope X Sn- The dependence of the obtained slope 
parameters on the charge number of the isotope is plot- 
ted in the inset (a) of Fig. [3] A second order polynomial 
function f{Z) was used to fit the results. Then we define 
an effective difference factor 



4// - cJ/fiz). 



(12) 



with f{Z) = 0.5848 + 0.04554Z - 5.356 x lO^^'Z^ being 
the fitted function. The reduced of the fit is 7.6 x 
10~®. From the dependence of deff on Sn as shown in 
the inset (b) of Fig. [3l def / always equals the neutron 
skin thickness for all the nuclei. If we could extract dgf / 
from experimental measurements, then the neutron skin 
thickness could be determined. To calculate de//, cf-N 
and (j'^i^ are necessary, crj^^ can be determined by the 
SAA model if we have the experimental proton radii data. 
Thus if (T-AT is measured experimentally, we could obtain 
deff- This indicates that Sn can be extracted from the 
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FIG. 2: (Color online). The N ~ Z dependence of the differ- 
ence factor d for Na, P, Ca and Ni isotopes with 5*^ varying 
from 0.0 fm to 0.5 fm. The horizontal lines just guide the eye. 
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FIG. 3: (Color online). The dependece of d with S„ for Na, P, 
Ca and Ni isotopes. The lines are the linear fit results. Inset 
(a) is the dependence of the slope of d with Z. The line is 
the fitted function f{Z). Inset (b) is the correlation between 
deff and Sn for Na, P Ca and Ni isotopes. The line is the 
function of y — x. For details see the text. 

experimental measurement of (T_7v- The accuracy of the 
extracted Sn is mainly determined by the uncertainty of 
CT-AT. Estimation based on Fig. [3] show that 5% error of 
a-N corresponds to around 0.05 fm error in Sn for Na 
isotopes. If the error of a-N could be smaller, then the 
uncertainty of the deduced Sn could become smaller also. 

Based on above discussion, a general phenomenon is 
shown that cr_Ar and 5'„ have strong linear correlation 
for neutron rich nuclei. Owing to the small uncertainty 
of Sn extracted from a-N data, a-N could be used as a 
observable to determine Sn- So we propose to use deff 
as a new observable to extract Sn experimentally. But 



it should be pointed out that Pygmy Dipole Resonance 
(PDR) exists in neutron rich nuclei, which has been in- 
terpreted as a resonant oscillation of neutron skin against 
the remaining isospin saturated neutron-proton core(25l- 
\2T\ . In this excitation mode, the Fermi distribution could 
not describe the nuclear densities of neutron rich nuclei. 
We are focusing on the relationship between Sn and a-N, 
similar correlation will also exist even if the shape of den- 
sity distribution is different since the RMS radii of neu- 
tron and proton are quantities describing the averaged 
distribution radii of the nucleon. Further investigations 
on the effect of the density shape over a-N and the un- 
certainty of the deduced Sn should be done. 

In summary, we have studied the relation between neu- 
tron removal cross section and neutron skin for neutron 
rich nuclei using SAA model. Strong linear correlation 
between a-N and the neutron skin thickness for neutron 
rich nuclei have been observed. Thus measuring a-N will 
be a very good method for studying the neutron skin 
size in nuclear surface. By defining an effective differ- 
ence factor, further analysis suggests that the neutron 
skin thickness could be extracted from the experimental 
measurement of a-N- We find that the uncertainty of Sn 
deduced from a-N data is quite small. So we propose to 
use def / as a new quantitative measure for neutron skin 
thickness in neutron rich nuclei. 
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